Delayed optical nonlinearity of thin metal films 
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Metals typically have very large nonlinear susceptibilities, whose origin is mainly of thermal 
character. We model the cubic nonlinearity of thin metal films by means of a delayed response 
derived ab initio from an improved version of the classic two temperature model. We validate our 
model by comparison with ultrafast pump-probe experiments on gold films. 



Noble metal nanostructures, such as thin films, grat- 
ings, multilayers and nanoparticles, have been largely ex- 
ploited in optics, especially within the fields of plasmonics 
and metamaterials. The linear optical response of these 
structures has been extensively studied both theoretically 
and experimentally, leading to the demonstration of un- 
precedented possibilities for extreme light concentration 
and manipulation (see [1| and references therein). Also, 
the intense nonlinear optical response of metallic nano- 
structures after intense excitation with fs-laser pulses has 
attracted increasing attention in view of the potential to 
achieve ultra- fast all-optical control of light beams 043 ■ 

Usually, the nonlinear response of the metal is mod- 
eled as a pure Kerr effect |2|-|4j], where the nonlinear po- 
larization is proportional to the cube of the electric field: 
Pnl = EoX^E 3 - in fact, this model though perfectly 
describing non resonant nonlincarities of electronic type 
that respond extremely fast (on the sub femtosecond time 
scale) to a driving electric field, turned out to be unsuit- 
able for fs and ps optical pulses. Actually, the values of 
the coefficient, measured by the z-scan technique, 
that can be found in the literature, differ by up to two 
orders of magnitude 0, [fij-Sl , clearly demonstrating that 
an instantaneous Kerr model is inadequate to describe 
the nonlinear response of metallic nanostructures. 

Pump-probe experiments in thin films (§-[IH and 
nanoparticles [l2|, reveal that the nonlinearity of met- 
als is due to the smearing of the electron distribution in- 
duced by intense optical absorption, resulting in a mod- 
ulation of the inter-band and intra-band transition prob- 
abilities with subsequent variation of the dielectric per- 
mittivity. The temporal dynamics of the system, which 
has been accurately interpreted according to the two- 
temperature model (TTM) [9|,|l4(, indicates that the non- 
linear response is dominated by a delay mechanism, but 
a theoretical formulation in terms of a non-instantaneous 
X^ 3 ) polarization is still lacking. 

In this Letter we derive, from an improved version of 
the TTM, a delayed third order non linear response suit- 
able for the description of optically thin metallic struc- 
tures. The outcome of our model is also quantitatively 
compared with experimental results from pump-probe 



spectroscopy on thin gold films. 

Starting from Maxwell equations (written in MKS 
units), neglecting transverse dimensions (i.e considering 
the propagation of plane waves), we can obtain the ID 
wave equation for the electric field E(z,t): 
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where c is the vacuum velocity of light, Sq is the vac- 
uum dielectric permittivity, e(w) = 1 + x( w ) an d x( w ) is 
the linear electric susceptibility (hat standing for Fourier 
transform). In the perturbative regime, the nonlinear po- 
larization can be expanded in Volterra series, accounting 
for small and non-istantaneous nonlinearity (lo| . Con- 
sidering only third order nonlinearity, nonresonant, in- 
coherent (intensity-dependent) nonlinear effects can be 
included by assuming the following functional form for 
the third-order polarization: 



PNL(z,t) = e 
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With the aim of quantifying the nonlinear response 
function by comparison with pump-probe experiments, 
we assume that the electric field is the sum of a powerful 
pump A and a weak probe B: E(z,t) — ^A(z,t)e lulat + 
~B(z,t)e iUht +c.c, with \A\ » \B\. In the slowly vary- 
ing envelope approximation (\d t A\ « \A\, \d t B\ « 
\B\), the evolution equation for the probe becomes: 
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where k\ = uj 2 Eb/c 2 = uj 2 e(ujb)/c 2 . In optically thin 
metallic structures we can assume that the pump is 
non-depleted, so that the nonlinear response is space- 
independent. By doing so, Eq. ([3]) can be written as 
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dz 2 



+ ^-[e b + Ae(\A(t-t d )\ 2 ,uj b )]B = 0, (4) 
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where the time-dependent nonlinear dielectric constant 
change Ae is a convolution between the pump pulse in- 
tensity and the third-order nonlinear response of the sys- 
tem (td being the pump-probe delay). Equation Q is 
a linear wave equation for the probe B, where the time 
delay t — td enters as a parameter. In the usual case of 
a short probe we can calculate the differential transmis- 
sivity (reflectivity) AT/T (AR/R) by considering t = 0, 
where the probe is peaked. 

The nonlinear response function X (*) can be derived 
from the extended two temperature model, that describes 
the evolution of the electrons and lattice temperature of 
a metal after absorption of a laser pulse Q : 



BT 

C e —^ = -g(T e -T l ) + aN 
at 

Q-^=g(T e -T l )+bN 



(5) 
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where C e and C; are the electronic and lattice heat ca- 
pacities, Ti is the lattice temperature, g is the electron- 
phonon coupling constant, N stands for the energy den- 
sity stored in the nonthermalized part of the electronic 
distribution, a is the electron gas heating rate and b is 
the electron-phonon coupling rate. P is the absorbed en- 
ergy density and is related to field intensity I through 
P{z,t) = {l-R- T)ae- az I{t); for a thin film, we can 
neglect the spatial dependence and assume a mean ab- 
sorbed energy density P{t) = (1 - R - T)I(t)/d (d, R, T, 
are the film thickness, reflection and transmission). The 
energy density N can be calculated as the convolution 
between the pump energy density and the thermaliza- 
tion response function hth(t) = exp[—t / T t h\H {t) (H(i) 
being the Heaviside function), with r t h = l/(a + b). The 
values a, b can be derived from a more sophisticated ex- 
tension of TTM proposed by Carpene [14(. It is worth 
noting that the two approaches give rather similar re- 



sults, but the rate equation approach Q is simpler. We 
obtain a = 1/2ti, t\ — E^tq / (Tiuji,) 2 and b — l/r ep , 
where Ep is the Fermi energy, tq — 128^/Stt 2 u> p (to p 
the plasma frequency), and r ep is the electron-phonon 
energy relaxation time. The value of To has been calcu- 
lated using the Lindhard dielectric function in the frame- 
work of Fermi liquid theory under the random phase ap- 
proximation [16]. This value turns out to be underesti- 
mated, mainly because of the d-band screening; a more 
correct value can be extracted from experimental data 

0. If AT e = T e - T a « T a , it is possible to con- 
sider the thermal capacity of electrons C e = 7T e = 7X0 
nearly constant: in this way the system (J5J) is linear 
and can be solved exactly. T can be estimated as a 
mean value between the minimum electronic tempera- 
ture T a and the maximum value, that can be estimated 

by T e _ max w *Jt% + 2/7 J P(t)dt. For most metals, in- 
cluding Au and Ag, C/ >> C e , allowing us to write a 
simple expression for AT e : 

AT e = / N(t- h) % + ( pf r e T' T Hit,)dh (6) 

J O e (pe + Cl )T t h 

1. e. the convolution between the energy density N 
and the two temperature system response h(t) — 

% + (C t +Ci)TtJ is now possible to calculate the 

total response of the extended two temperature system 
as the convolution of h t h and h: 
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This way the electronic temperature change is simply cal- 
culated as the convolution between the pump and the to- 
tal response: AT e = J P(t — ti)htot(ti)dtx- Assuming a 
gaussian pump P(t) — Pq exp[— 2(t/T ) 2 } we obtain 



AT„ 
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Following the same procedure a similar expression can be 
obtained also for the lattice temperature variation ATl 
and the non-thermalized electron energy density N. 

We tested the validity of our approach by studying the 
temperature changes of a thin gold film (d = 20 nm) 
deposited on a sapphire substrate (n s = 1.75), after irra- 
diation with a pump pulse of A a = 950 nm wavelength, 
T = 140 fs duration and F — 30 ^J/cm 2 fluence. Pa- 



rameters for the TTM were taken from literature [ll|, [14 1 
as T a — 300 K, 7 = 70 J i 



k-2, a = 2.5 • 10 6 j 
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K" 1 , g = 2- 10 16 W m~ 3 K^ 1 , r ep = 1.4 ps, r = 6.5 fs 
16], Ep — 7.3 eV. Figure [T] shows the results of numeri- 
cal solution of Eqs. (0, of numerical solution of Carpene 
model [14j and of analytical formula (j8]). The agreement 
between analytical approximation and TTM is remark- 
able, indicating that in this range of electronic temper- 
ature the assumption of a constant heat capacity is rea- 
sonable. Carpene model and TTM give practically the 
same outcome for what concerns lattice temperature, but 
a small discrepancy can be noticed for what concerns the 
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FIG. 1. (color online) Temperature dynamics after irradi- 
ation of a thin gold film deposited on a sapphire substrate, 
from TTM of Eq. ((5| (dashed black curves), from Carpene 
model [14fl (dotted blue curves), and from analytical formula 
of Eq. ([8]) (solid red curve). Inset: zoom-in of the initial tem- 
perature dynamics; normalized energy density stored in non 
thermalized electrons (dash-dot black curve); pump intensity 
profile (thin dashed green curve). See text for parameters. 



peak of electronic temperature. This feature is mainly 
due to the quite long thermalization time T t u and will 
produce a slightly higher estimation of the maximum di- 
electric constant changes. In this figure we also show the 
time evolution of the energy density stored in the non 
thermalized electrons N. It acts like a delayed effective 
pump for electron and lattice, slowing the rise of AT e 
and AT;. The dynamics of N is faster than the ther- 
malized one and it is responsible for the small and quick 
dielectric permittivity changes observed when probing in 
a frequency range far from inter-band transitions [9(. 

Variation in electronic temperature induces variation 
in the dielectric constant through inter-band (bound elec- 
trons) transitions Q. Additional contributions to the 
change in the dielectric constant are given by intra-band 
transitions (that depend on lattice temperature) and by 
non thermalized electron energy distribution. These ad- 
ditional contributions are much smaller than the former 
and can be resolved only in spectral regions where inter- 
band transitions are not efficient [9|. In order to obtain 
a reasonably simple model that however correctly grasps 
all the relevant phenomena, from now on we concentrate 
on inter-band transitions of thermalized electrons. 

Energy injection into the conduction band smears 
the electron distribution around Ep with reduction (in- 
crease) of the occupation probability of the electron 
states below (above) Ep. As a consequence, a modula- 
tion of inter-band transition probability is induced, with 
increased (decreased) absorption for transitions involv- 
ing final states below (above) Ep. Inter-band transitions 
in noble metals are dominated by d-band to conduction 
band transitions near L and X points in the irreducible 
zone of the Brillouin cell. In the constant matrix element 
approximation, the variation of the imaginary part of the 
inter-band dielectric function of gold can be computed as 



Ae 2 (to b ,T e ) 
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where m is the free electron mass, Pl and Px are the 
electric-dipole matrix elements, and AJl and AJx the 
temperature induced variations of the Joint Density of 
States (JDOS) for o?-band to conduction band transitions 
near L and X respectively. Such variations can be com- 
puted as fl7j ]: 



AJ L>x (ub,T e ) 



D L , x {E,uj b )AP{E,T e )dE (10) 



where Dlx{E, uib) is the Energy Distribution of the Joint 
Density of States (EDJDOS) of the considered transi- 
tions (with respect to the energy of final state E), and 
AP(E, T e ) = P{E, T e ) - P(E, T a ) is the temperature in- 
duced variation of the occupation probability for the final 
state (being P(E, T e ) = 1 - f(E, T e ) the probability that 
the final state is empty and f(E, T e ) the Fermi-Dirac 
function). For small temperature changes AT e , the fol- 
lowing approximation holds: 
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FIG. 2. (color online) Spectral variation of the (a) imaginary 
part and (b) real part of the inter-band dielectric function of 
gold for three different values of AT e , computed from Eqs. ([!)]■ 
Illl) . Linearized temperature dependence provided by Eq. pip 
(solid lines) is compared with exact solution (dashed lines). 

The EDJDOS was numerically computed under 
parabolic band approximation following the approach 
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described in [FjJ, taking effective masses, energy gaps, 
dipoie matrix eiements and integration limits E' L x and 
E'l x as reported in [l8[. The Ae 2 resulting from Eqs. ([9} 
as a function of photon energy for three different tem- 
perature variations is shown in Fig. 2(a). The variation 
of the real part of the inter-band dielectric function Aei 
computed by Kramers-Kronig analysis of Ae 2 is reported 
in Fig. 2(b). Note that the approximation provided by 
Eq. (fTTj) is accurate for AT e as high as 100 K. 

For a quantitative validation of our model, we com- 
puted the transient differential reflection (transmission) 
AR/R (AT/T) of the 20- nm thin gold film according to 
standard thin film formulas [l9[ and Eqs. (|8l fTT|) . Results 
are reported in Fig. 3, compared with experimental data 
taken from the literature [9j and full-numerical computa- 
tion from the extended TTM. The probe photon energy 
is chosen at the maximum response of the thin film [9j, 
i.e. huj b = 2.48 eV for AR/R and hu) b = 2.43 eV for 
AT/T. We found that the analytical model is in quanti- 
tative agreement with experimental data, with only slight 
deviations in the long-time range (t > 4 ps) where the 
contribution from the lattice is dominant. 




1 2 3 
Probe delay (ps) 

FIG. 3. (color online) Transient (a) reflectivity and (b) trans- 
missivity vs probe time delay, accounting for pump-probe 
cross-correlation (shown as dashed green line). 



According to above analysis, a delayed cubic response 
of the metallic medium can be written as follows: 



X (3) (t) = ^acne h tot (t) 



(12) 



The magnitude of the nonlinear response can be esti- 
mated by assuming continuous-wave pump and probe. 
In this case the nonlinear response is much faster than 
the fields, and the convolution reduces to the integral of 



X^[t). Assuming C; >> C e we obtain: 

-m r \ dAe(uj b ) ar th T 
X K >{u} b ;u} a ,-uJa,Ub) = — kf^ a(uj a )cn(uj a )e - 



(13) 

(the hat is added because assuming c.w. fields is equiv- 
alent to give a frequency domain description of X )■ 
As an example, in the case of an infrared pump A a = 
950 nm and visible probe A& = 500 nm, we obtain 
X (3) (w fc ;w a ,-w a ,w fc ) « (-8.4+ Hi) • 10- 8 esu= (-1.2 + 
1.5i) ■ 10~ 15 m 2 /V 2 . This huge value (six order of magni- 
tudes greater than in fused silica) is due to the resonance 
of the probe with interband transitions of gold. There- 
fore, it should be stressed that this kind of nonlinearity is 
strongly dispersive, and can be assumed constant only for 
reasonably band-limited pulses centered around Lo a and 
u)]). Changing either pump or probe frequency results in 
a drastic change of x^H^b), confirming that this kind of 
third order nonlinearity is not of Kerr-type. 

To conclude, we introduced a theoretical model for the 
delayed nonlinear response in optically thin noble-metal 
structures. A non-instantaneous x^ coefficient is de- 
rived from an extended version of the TTM and semi- 
classical theory of optical transitions in solids. Our theo- 
retical predictions turned out to be in quantitative agree- 
ment with experimental results from pump-probe exper- 
iments in thin gold films. 
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